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Abstract. We will introduce a new method to derive the orbital elements of the 



spectroscopic binary stars. Fitting the adopted radial velocity curves on the mea- 
sured experimental data, the corresponding radial acceleration data are obtained. 
The orbital parameters are estimated by the nonlinear least squares of the radial 
velocity-acceleration curves. This method is implemented for the three double line 
spectroscopic binary systems RZ CAS, CC Cas and VI 130 Tau. Numerical cal- 
culations show that the result are in good agreement with the those obtained by 
others using the method of Lehmann-Filhes. 
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1. Introduction 

Determining the orbital elements of the binary stars help us to obtain the necessary 
information such as the mass and the radius of stars which play the important roles during 
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the evolutions of the stellar structures. Analyzing both the light and the radial velocity 
curves deducing from the photometric and the spectroscopic observations, respectively, 
yields to derive the orbital parameters. One of the usual methods to analyze the velocity 
curve is the method of Lehmann-Filhes (Smart 1990). Here we introduce a new method 
to derive these parameters by the nonlinear regression of the radial velocity-acceleration 
curves. We test our method for the three double-line spectroscopic binary systems RZ 
Cas, CC Cas and VI 130 Tau which have the following properties. 

The Algol system RZ Cas is the semi-detached eclipsing binary stars with a period 
P = 1.195 days. There is a mass transfer between the two components. Some evidences 
existed for periodicity in the period changes due to apsidal motion. RZ Cas is also a radio 
source and an X-ray source. See Riazi et al. (1994) and Maxted et al. (1994) and references 
therein. CC Cas is the massive O-type binary system with a period P — 3.366 days. This 
system is a detached binary with the two components within the main sequence band. 
It is an actively interacting system. The system shows a mass loss rate 3 x 10 m0yr _1 
due to the stellar winds. Such mass loss from the system yields to a change in the orbital 
period by the rate of less than 3 x 10 -6 days per century (Hill et al. 1994). V1130 Tau 
is a very closed detached system including the two components with similar size and 
brightness. It consists of two almost identical F3 V stars on a tight orbit with a short 
period P = 0.799 days. (Rucinski et al. 2003). 

In Sect. 2, we reduce the problem to solving a equation for the radial velocity- 
acceleration relation which is a nonlinear function in terms of the orbital parameters. In 
Sect. 3 the nonlinear regression technic for estimating the orbital elements is discussed. 
In Sect. 4 the numerical results implemented for the three different binary systems are 
reported. Section 5 is devoted to concluding remarks. 

2. Formulation of the problem 

The radial velocity of star in a binary system is defined as follows 

V r = V cm + Z, (1) 

where V cm is the radial velocity of the center of mass of system with respect to the sun 
and 

Z = K[cos(9 + w) + ecoscj], (2) 

is the radial velocity of star with reference to the center of mass of the binary (see Smart 
1990). Note that the dot denotes the time derivative. In Eq. (J2J, 6, lo and e are the angular 
polar coordinate, the longitude of the perihelion and the eccentricity, respectively. Also 
2ir asini 

where P is the period of motion and the inclination i is the angle between the line of 
sight and the normal of the orbital plane. 
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The radial acceleration of star can be obtained as 



Z = -Ksm(0 + lo)9. (4) 

Using Kepler's second law and the relations obtaining for the orbital parameters in the 
inverse-square field as 

6 = h/r 2 , (5) 
2tt 



h= T a 2 Vl^, (6) 

r= (7) 
l + ecos0' 1 1 

the radial acceleration, Eq. (0J, yields to 
—2-kK 

Z = p (1 _ e 2)3/2 M<> + + 6COS0) 2 , (8) 

where r, a and h are the radial polar coordinate, the semi major axis of the orbit and 
the angular momentum per unit of mass, respectively. 

Substituting for 6 from Eq. in Eq. (J8J, the result reduces to 

PZ = ( 1 _ e 2 )3 /2 sin (cos (Z/Jf-ecosw)) 

x|l + ecos ( - w + cos _1 (Z/X - ecosw))| . (9) 

To simplify the notation further, we let Y = PZ and X = Z. Hence Eq. JJJJ describes a 
nonlinear relation, Y" = Y(X, K, e,u>), in terms of the parameters K, e and u>. 

One can show that the adopted spectroscopic elements are related to the orbital 
parameters. First according to the definition of the center of mass, the mass ratio in the 
system is obtained as 

(10) 



m s dp sin i 

where subscripts p and s return to the primary and the secondary components of the 
system. 

Then from Kepler's third law and Eq. 1)10(1 . the result reduces to 

. 3 . . .,a p smi + a s smi 2 
m p sin j = a 8 sint(— — j , (11) 

where a, P and m are expressed in AU, years and solar mass, respectively. The similar 
relation is obtained for the secondary component only by replacing p to s and vise versa, 
in Eq. (|1 lfl . Note that in Eqs. (110ft and 1)11(1 the parameter asini is related to the orbital 
parameters by the aid of Eq. (j2J. 
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3. Nonlinear regression of the radial velocity-acceleration curve 

To obtain the orbital parameters K, e and uj in Eq. @, we use the nonlinear regression 
method. In this approach, the sum of squares of errors (SSE) for the number of N 
measured data is calculated as 

N N 

SSE = Y J {Y*-Y l ) 2 =^[F 4 -r(X l; X, e ,c)] 2 , (12) 

z=l i=l 

where Yi and Yi are the real and the predicted values, respectively. To obtain the model 

parameters, the SSE should be minimized in terms of K, e and u> as 

dSSE dSSE dSSE . . 

~3K - = — = — =°- < 13 > 

To solve Eq. JT^jl , we use the SAS (Statistical Analysis System) software. Note that the 
nonlinear models are more difficult to specify and estimate than linear models. For in- 
stance, in contrast to the linear regression, the nonlinear models is very sensitive to the 
initial guesses for the parameters. Because in practice, SSE may has to be minimized in 
several points in the three dimensional parametric space including K, e and to. However 
the final goal is finding the absolute minimum. Hence choosing the relevant initial pa- 
rameters yields to the absolute minimum of SSE which is also stationary. That means 
if one change the initial guesses slightly, then the result reduces to the previous values 
for the parameters. But if the regression converged at the local minimum, the model 
would not be stationary. See Sen & Srivastava (1990) and Christensen (1996). To avoid 
the problem, two things can help us, first by taking into account the bound intervals for 
the parameters which causes the model to converge quickly. For instance one can put 
(0 < e < 1) for the binary stars which have closed orbits. The last one is checking the 
differences between the observational data and the fitted diagrams using Eq. (JHJ. 



4. Numerical results 

Here we test our new method to derive the orbital and the combined elements of the 
binary stars. To do this we consider the three different double line spectroscopic systems, 
RZ cas, CC Cas and VI 130 Tau. Using the measured experimental data for the radial 
velocities of two components of these systems obtained by Maxted et al. (1994) for RZ 
Cas, Hill et al. (1994) for CC Cas and Rucinski et al. (2003) for VI 130 Tau, the suitable 
fitted velocity curves are plotted in terms of the photometric phase in Figs. ^ 01 an d 
[7| The best fits are obtained by the polynomials of the intermediate degrees, e.g., sixth 
and seventh degrees. Because the low degrees do not suitably cover the data and high 
degrees take the fluctuations of the data which yields to the divergency in the nonlinear 
regression of Eq. ©. The velocity of the center of mass of the stars are obtained by 
calculating the areas up and down of the radial velocity curves. Everywhere these areas 
become equal to each other then the velocity of the center of mass is obtained. 
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The radial acceleration data corresponding with the radial velocity data are obtained 
by taking the derivative of the adopted radial velocity curves. Figures [21 El EI El El and |^1 
show the radial acceleration scaled by the period versus the radial velocity data for the 
primary and the secondary components of RZ cas, CC Cas and VI 130 Tau, respectively. 
The solid closed curves are the result of the nonlinear regression of Eq. , which their 
good coincidence with the observational data yields to derive the optimized parameters 
K, e and lj. 

The orbital parameters obtaining from nonlinear least squares of Eq. (|12|1 for RZ 
Cas, CC Cas and VI 130 Tau are tabulated in Tables ^ El El respectively. Table E] shows 
that for RZ Cas: 1) V cm and K p for the primary component are in good agreement with 
the those obtained by Maxted et al. (1994) and Duerbeck & Hanel (1979). 2) e and cu 
are closed to the results obtained by Horak (1952). 3) V cm and K s for the secondary 
component are coincide with the results derived by Maxted et al. (1994). Whereas for e 
and w, we had no any references for comparing. 

TableEJshows that for CC Cas: 1) V cm and K for the primary component are in good 
analogy with the results obtained by Hill et al. (1994). For the secondary component, K 
is in relatively well coincidence with the those derived by Pearce (1927) and Hill et al. 
(1994). But for V cm is not, however its interval is compatible with the result of Hill et al. 
(1994). The difference may be caused from the existence of high scattering for the data 
around the maximum and minimum of the radial velocity curve (see again Fig. EJ) ■ 

TableEJshows that for V1130 Tau, 1) for both components, K is in very good concord 
with the result estimated by Ruciniski et al. (2003). 2) For V cm , its average value for the 
two components is in good agreement with the results of Rucinski et al. (2003). 

The combined spectroscopic elements including m p sin 3 i, m s sin 3 i 7 (a p + a s ) s'mi and 
m p /m s are calculated by substituting the estimated parameters e, K and u> in Eqs. Ipj|. 
H10(l and JUJ. The results obtained for the same three previous systems are tabulated 
in Tables El El El Tables show that our results are in good agreement with the those 
obtained by Maxted et al. (1994), Hill et al. (1994) and Rucinski et al. (2003) for RZ 
Cas, CC Cas and VI 130 Tau, respectively. 

5. Concluding remarks 

A new method to derive the orbital elements of the spectroscopic binary stars is in- 
troduced. These parameters are obtained from the nonlinear regression of the radial 
velocity-acceleration curve of the system. Using the measured experimental data for ra- 
dial velocities of RZ Cas, CC Cas and V1130 Tau obtained by Maxted et al. (1994), Hill 
et al. (1994) and Rucinski et al. (2003), we find the orbital elements of these systems 
by the mentioned method. Our numerical results show that the results obtained for the 
orbital elements and the combined spectroscopic parameters are in good agreement with 
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the those obtained by others via the method of Lehmami-Filhcs. In a subsequent paper 
we intend to test numerically our method for the other different systems. 
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Fig. 1. The radial velocities of the primary and secondary components of RZ Cas plotted 
against photometric phase. The observational data have been measured by Maxted et al. 
(1994). Legend is given in the right hand up caption. 
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RZ-cas (primary) 

1500 i 1 1 1 1 



1000 ■ 
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Fig. 2. The radial acceleration scaled by the period versus the radial velocity of the 
primary component of RZ Cas. The solid curve is obtained from the nonlinear regression 
of Eq. © . The plus points are the experimental data. 
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RZ-cas (secondary) 
1500 i 1 1 1 1 1 1— 
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Fig. 3. Same as Fig. for the secondray component of RZ Cas. 
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CC-cas 
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Fig. 4. Same as Fig. Q f° r CC Cas. The observational data are belong to Hill et al. 
(1994). 
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CC-cas (primary) 
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Fig. 5. Same as Fig. for the primary component of CC Cas. 
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CC-cas (secondary) 
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Fig. 6. Same as Fig. for the secondary component of CC Cas. 
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V1130-Tau 
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Fig. 7. Same as Fig.Q for VI 130 Tau. The observational data have been deduced from 
Rucinski et al. (2003). 
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V1130-Tau (primary) 
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Fig. 8. Same as Fig. for the primary component of VI 130 Tau. 
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V1130-Tau (secondary) 
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Fig. 9. Same as Fig. for the secondary component of V1130 Tau. 
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Table 1. Spectroscopic orbit of RZ Cas. Values bold type are assumed. 





This Paper 


Maxted et al. (1994) 


Hork (1952) 


Duerbeck & Hanel (1979) 


Primary 










V cm (ferns -1 ) 


-45.08 ± 0.93 


-45.5 ±0.7 


40 ±2 


-45.5 ± 1.1 


Kp (fcms -1 ) 


70.4 ± 0.34 


70.9 ± 0.8 


68 ±2 


70.5 ± 1.4 


e 


0.013 ±0.002 


0.000 


0.01 ±0.03 


0.024 ±0.023 


o 


250 ± 11.5 


— 


240 ±160 


87 ±2 


Secondary 










Vcm (fcms -1 ) 


-40.68 ± 3.69 


-41 ±4 








213.35 ±0.45 


213 ±4 






e 


0.008 ±0.001 


0.000 






o 


88.84 ± 7 









Table 2. Combined spectroscopic orbit of RZ Cas. 



Parameter This paper Maxted et al. (1994) 



m p sin 3 i/Mq 


2.126 ±0.002 


2.16 ±0.07 


m s sin 3 i/Mq 


0.701 ± 0.005 


0.715 ±0.02 


(a p ± a s ) sin i/10 6 fcm 


4.663 ± 0.007 


4.68 ±0.07 


nip/ms 


3.03 ±0.02 


3.02 ±0.07 



Table 3. Same as Table Q] for CC Cas. 



This Paper Hill et al. (1994) Pearce (1927) 

Primary 

V cm (fcms -1 ) -7.23 ±2.68 -7.3 ± 4.2 -9.9 ± 2.4 

Kp{kms~ 1 ) 117.87 ±0.51 116 ±1.8 138 ± 4 

e 0.0025 ± 0.0005 

lo(°) 330 ±2 
Secondary 

V cm (kms- 1 ) -8.63 ±7.31 -12.3 ±6 3.2 ± 7 

K s (ferns' 1 ) 285.63 ± 1.14 279.4 ± 2.3 288 ± 9 

e 0.023 ±0.005 

lo{°) 54 ±13 



Table 4. Same as Table El for CC Cas. 

Parameter This paper Hill et al. (1994) Pearce (1927) 

m p sin 3 i/M e 16.2 ± 0.2 15.3 ± 0.3 18.3 ±1.4 

m s sin 3 i/M 6.69 ±0.08 6.3 ± 0.2 8.8 ± 0.6 

(a p + a a )smi/R Q 26.83 ± 0.08 26.2 ± 0.2 28.3 ± 0.6 

m p /m s 2.42 ±0.02 2.41 ± 0.04 2.09 ± 0.09 
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Table 5. Same as Tabled f° r VI 130 Tau. Number in parenthesis are errors in the final 
digits. 





This Paper 


Rucinski et al. (2003) 


Primary 






V / 


-14.45 ±0.89 


-12.74(0.46) 


K v \ kms~] 


146.61 ± 1.01 


147.21(0.63) 


e 


0.005 ± 0.003 




w(°) 


126 ± 42 




Secondary 






Vcm (fcms" 1 ) 


-10.9 ±0.9 


-12.74(0.46) 




159.36 ±0.61 


160.11(0.74) 


e 


0.004 ± 0.004 






292 ±7 





Table 6. Same as Tabled for V1130 Tau. 



Parameter 


This paper 


Rucinski et al. (2003) 


m p sin 3 i/Mq 


1.23 ±0.02 




m s sin 3 i/Mq 


1.13 ±0.2 




(cip + a a ) smi/R Q 


4.82 ±0.02 




m p /m s 


0.919 ±0.009 


0.919(7) 


(m p + m s ) sin 3 i/Mq 


2.37 ±0.02 


2.408(32) 



